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ABSTRACT 

Electron temperatures measured by electrostatic probes and radar 
backscatter are distinct physical quantities with the temperature from 
each technique determined from a different moment of the electron 
distribution function. Numerical inequality of temperatures results 
from a non- Maxwellian electron distribution function or equivalently, a 
non-equilibrium electron plasma. Probe and backscatter electron tem- 
peratures are studied for low and high energy (isotropic) distortions of 
the distribution function. The non-equilibrium plasma generally pro- 
duces higher probe than backscatter temperatures, however the tem- 
perature difference is small for distortions due to realistic photoelectron 
fluxes. If large temperature differences occur in the ionosphere, both 
probe and backscatter temperatures would differ from the temperature 
characterizing the average electron kinetic energy, and a single tem- 
perature applicable to a variety of physical processes would no longer 
exist. 
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INTRODUCTION 


Comparisons of electron temperatures reduced from electrostatic probe 
and radar backscatter measurements have resulted in a temperature discrep- 
ancy, with higher probe than backscatter temperatures [Hanson et al, 1969; 
Carlson and Sayers, 1970; Booker and Smith, 1970; Brace and McClure, 1971] . 
Some comparisons have demonstrated agreement among the probe and back- 
scatter temperatures [Brace, et al, 1969; Taylor and Wrenn, 1970], Two possible 
explanations for the discrepancy have been suggested recently: the presence of 
suprathermal electrons in the electron distribution function (Waldteufel, 1969) 
and the slow variation of the contact potential of the probe surface (Carlson and 
Sayers, 1970). The contact potential variation in laboratory probes has been 
eliminated by probe heating and by a rapid sweep of the probe potential 
[Waymouth, 1959], However it has not been firmly established that rocket and 
satellite probe electron temperatures are significantly affected by the contact 
potential variation. Several authors have discussed this and other possible 
sources of error in the probe temperature measurements and have in general 
confirmed the technique [Wilmore, 1970; Brace et al. , 1970]. 

As regards the effect on the probe and radar temperatures due to an en- 
hanced energetic electron population, only the effect on the radar temperatures 
has been discussed in the literature. 



t 


We investigate the effect of an isotropic non-Maxwellian electron distribution 
function on both probe and radar electron temperatures. This distribution repre- 


sents a simplification of actual ionospheric conditions, nevertheless it allows us 
to study the important case of energetic photoelectrons as well as to illustrate the 
basic theoretical difference in the temperature parameters of the two techniques. 
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DISTRIBUTION FUNCTION 


The isotropic electron distribution, f(v), is used in two representations: 

1) as a sum of two Maxwell distributions at temperatures Tj and T 2 and with 
populations 1-p and p respectively, 


f (v) = (i _ p) 
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and 2) as a Laguerre polynomial expansion 
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The two-temperature distribution, Eq. 1, may be used to represent the presence 
of energetic photoelectrons when T 2 > T , and P < < 1, and a distorted low energy 
electron population for nearly equal values of Tj and T 2 . The coefficients A2 , 

A3, etc. in the polynomial representation determine the deviation from a Maxwell 
distribution. WhenT 2 /Tj < 10, then the two temperature distribution can be 
represented by the Laguerre polynomial expansion, however the latter can repre- 
sent a greater variety of non- equilibrium situations. 

TEMPERATURE DEFINITIONS 

The probe and radar temperature parameters are determined from different 
moments of the electron distribution function. Each temperature is therefore a 
distinct physical parameter since it samples a different portion of the electron 
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energy spectrum. There is strict numerical equality of the temperatures only 
when the distribution function is pure Maxwellian. Historically the probe and 
radar temperatures have been considered identical physical parameters since 
the assumption of a Maxwell distribution has been inherent in the theories of the 
two techniques. We shall examine the extent of the probe and radar temperature 
difference for an assumed deviation from a Maxwell distribution, or degree of 
departure of the ionosphere from equilibrium. In the following we review the 
definitions of four instrument temperatures (three probes and one radar) and 
the thermodynamic temperature, which definitions are valid for any isotropic 
distribution function (with unit normalization). 

The thermodynamic temperature, T th , is determined from the average 
energy of the entire electron population, 

3/2 k T th = Jd V 1 m V 2 f (v). ( 3 ) 

The radar backscatter electron temperature, T b , is one of several param- 
eters determining the shape of the backscatter power spectrum [Rosenbluth 
and Rostoker, 1962 ] . The ion component of the power spectrum, neglecting 
terms of order (m e /m.) 1/2 , depends on the electron distribution function only 
through the integral, 


kT b 



(4) 
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Thus the effect of a distorted isotropic electron distribution function on the radar 
electron temperature can be studied without analyzing the power spectrum or 
the correlation function. The velocity moment defining T gives more weight 

b 

to the low energy electrons than the moment defining T th . 

The three probe temperatures are defined in terms of the retarded electron 
current 


i = K 



v d v 



2 e v\ 
m / 


f (v). 


( 5 ) 


where K is a constant depending on the probe geometry and V is the absolute 
magnitude of potential difference between the probe and plasma. The probe 
temperature, T plate , corresponds to the temperature measured by the a.c. mode 
Langmuir plate [Wrenn, 1969] . This device, by operating in the a.c. mode, 
measures the ratio of the first to the second derivative of the electron current. 
Thus T plate is defined by, 


T 

p 


late 



vdv f (v)/f 



( 6 ) 


The extent of the variation of T , afe with potential can in principle provide in- 
formation on the electron distribution function, even though it does not directly 
measure f(v). 

The temperature, T eyl , corresponds to the temperature measured by the 
cylindrical electrostatic probe [Findlay and Brace, 1969] . The temperature 


4 


T cyl is that parameter which gives the best fit of the retarded electron current 
to an exponential of the potential over a given range of potentials. In order to 
separate the retarded current from the net current, we adopt the procedure set 
forth by Theis [1968] . The parameter T cy] is determined from a best fit to 
the net current, using the four parameter expression, 


net 


= C, e 


-eV/kT 


eyl 


c 2 ♦ c 3 V. 


(7) 


Some information on the electron distribution function can be obtained by varying 
the voltage interval over which the fit is made. 

The temperature, T trap , corresponds to the temperature measured by the 
planar electron trap [Donley, 1969] . This device measures the retarded 
electron current over a wide voltage range, enabling it to detect both the thermal 
and the energetic electrons. We define T tr . p from the slope of In i versus V, 




v d v f (v) 


( 8 ) 


2eV/m 


We now specialize the five temperature definitions for the case of the two- 
temperature distribution function. The thermodynamic and backscatter tem- 
peratures have a simple form, 


T th = (l-P) T. +P T 2 

(9) 

1/T b = (1 -P)/T, +P/T 2 . 

(10) 
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Asa consequence of these definitions, T 2 T fc with equality holding when p = 
0, 1, and with the greatest difference occuring when p ■ 1/2. The retarded 
electron current reduces to the sum of two exponentials, 


(1 -P) 



-eV/kT 
e 1 


+ P 


-eV/kT, 

e 


( 11 ) 


thus the trap and plate temperatures are given by, 

^tr«p = (l ~ q) Tj + q Tj, 

1/T pi.t. = (1 “< 0^1 +q/T a . 


( 12 ) 

(13) 
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It follows from Eqs. 12-14 that for a two- temperature distribution, T trap 2 . 

T p i ate , with equality holding for p = 0, 1. (This inequality holds only if both 
measurements are taken at the same voltage.) The parameter, T cyl is given by 
fitting the current of Eq. 7 to the form of Eq. 11. 


Next we specialize the temperature definitions for the Laguerre polynomial 
expanded distribution function. The thermodynamic and backscatter tempera- 
tures are, 

T = T. (15) 
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(16) 


T„ = T/ £**• 

-£■ o 

There is no inequality relating T th and T b and either temperature could be 
larger than the other. The retarded electron current, used to define the probe 
temperatures, is given as, 


/f e -eV/kT L A* 

'C* o 



( 17 ) 


The temperature T cyl is obtained from the four parameter fit to this current 
expression, and the trap and plate temperatures are 



We have illustrated the formal difference between the five temperature 
parameters, T th , T b , T pI>t# , T ey , , and T ; in following sections we shall 
present examples of numerical temperature differences for given distortions 
of the distribution function from a Maxwell distribution. And in order to connect 
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the distortion of the distribution function with the degree of departure of the 
ionospheric electrons from equilibrium we shall include an evaluation of the 
degree of departure parameter, L, [Lunc, 1963 ]. The parameter L is defined 
by, 


•I 


L + 1 = I d v f 2 /f eq, 


( 20 ) 


where f is the equilibrium distribution function. The state of the electrons is 
considered far from equilibrium when L > > 1, close to equilibrium when L < < 1, 
and in a transition regime when L ~ 1. 

The voltages used in computing the probe temperatures are: .5 volts for T 

plate 

.1 volts for T trap ; .5 to 2.5 volts for T . 

PHOTOELECTRONS 

The effect on the five temperature parameters of energetic photoelectrons 
is examined with the two-temperature distribution function. The superthermal 
electrons are represented by a Maxwell distribution at temperature T 2 and with 
relative density p, p « 1, T 2 > T , where T is the temperature of the low 
energy or thermal electrons. A convenient description of the photoelectrons is 
to specify their temperature, T 2 , and normalized flux, Photo, 

Photo = p /T^. (21) 


(The flux is given by, n i k/27rm Photo.) 


The range of values of T 2 and Photo which characterize the ionospheric 
photoelectrons are estimated from three sources. In the first, Dalgarno McElroy 
and Stewart computed values of the equilibrium photoelectron flux [ Dalgarno 
et al., 1969] . A Maxwellian fit to their 350 km flux values yields the values, T 2 
%70,000°K and Photo £ 10 4 /n, where n is the total electron number density in 
part. cm* 3 . The electron trap data of Donley at 1410 km yield the approximate 
values T 2 ~ 90,000°K and Photo ~.8 [Donley, 1969]. And finally Huang has 
made rocket measurements of superthermal electrons in the altitude range from 
120 km to 240 km [Huang, 1969] . The low altitude values are approximately, T 2 ~ 

10 ,000 °K to 30 ,000°K and Photo £ .3 to .5. 

The above values suggest that the ionospheric photoelectrons are charac- 
terized by the range of values, 

T 2 = 10,000 °K - 100.000 °K 
Photo = .1 - 1 

The three sources are displayed as shaded areas in Figure 1 bounded by these 
ranges. The differences among the temperature parameters are displayed in 
Figures 2-5 using the same coordinates as in Figure 1 so that they may be 
correlated with the photoelectron regions. 

The radar backscatter temperature, T b , differs from the thermal electron 
temperature, Tj , by no more than 1% over the region of Figure 1, therefore 
for photoelectrons 

T - T . 

*b V 
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In Figures 2-5 are displayed contours of the relative differences between 
the temperatures T th , T piata , T trap , and i respectively with the thermal 
population temperature T t = 1500°K. The relative differences decrease as Tj 
increases (or as the electrons approach equilibrium). Figure 2 shows that the 
temperature difference between T th and Tj is greater than 5% for Photo > .5, 
but that it does not exceed 20% over the entire range. The temperature difference 
between T plate and T, in Figure 3 does not exceed 3% for the three regions of 
Figure 1. The difference between T and T remains less than 5%. In some 

trap 1 

regions the inequality, T trap ^ T plate » is violated because the voltages at which 
the temperatures are defined were chosen to reflect the conditions of the experi- 
ments. The difference between T cyl and Tj in Figure 5 varies between 5% and 
10% in the low altitude region and is less than 3% in the higher altitude regions. 

Figure 6 illustrates the variation of the five temperatures with the photo- 
electron energy for a fixed flux, Photo = 1, and for T = 1500°K. 

We conclude that photoelectrons produce no significant temperature differences 
among the four instrument temperatures, T. , T , T , , and T . . An 

evaluation of the parameter L, Eq. 20, specifying the degree of departure of the 
electron gas from equilibrium shows that the electrons are far from equilibrium. 
(To obtain a finite value of L it is necessary to introduce a cutoff energy E above 
which there are no photoelectrons, a typical calculation yields L = 3 x io 6 for 
T t = 1500°K, T 2 = 30,000°K, Photo = .5, and E = 10 ev.) Thus the highly non- 
equilibrium state of the ionosphere as produced by photoelectrons, although it 
causes such phenomena as dayglow, yet it does not affect the probe and radar 
temperature measurements; they essentially measure only the thermal electron 
population. 


The neglect of anisotropy in the photoelectron velocity distribution would 
not alter this conclusion. A particle precipitation producing energetic electrons 
at fluxes substantially higher than the fluxes of Figure 1 would produce elevated 
probe and radar temperatures as well as temperature differences among the 
five temperature parameters. Large increases in the cylinder probe tempera- 
ture have been observed in coincidence with soft electron flux precipitating in 
the auroral region [Brace and Findlay, 1970]. 

DISTORTED THERMAL ELECTRONS 

The effect of a distorted thermal (low energy) electron population on the 
probe and radar temperature parameters is first evaluated using the two- 
temperature distribution function. We choose the values Tj = 1000, P = .5 (to 
maximize the distortion) and the range of values T 2 /T = 1 to 2. The resulting 
temperatures are displayed in Fig. 7. For the range of parameters shown, the 
three probe temperatures are greater than T th , while the radar temperature is 
smaller than T th . Therefore whenever there is a significant probe - radar 
temperature difference, then both instrument temperatures differ appreciably 
from the thermodynamic temperature. 

The measure of the degree of deviation from equilibrium, L, varies from 
L = .0031, .011, .023 at T/T, = 1.1, 1.2, 1.3 to L = .133, .153 at T 2 /T, = 1.9, 
2.0 respectively, thus we classify the two-temperature distribution used in Fig. 

7 as lying in the transition regime or moderately close to equilibrium. This 
demonstrates that moderate to large temperature differences do not necessarily 
require the ionosphere electrons to be far from equilibrium. 
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Finally we examine the effect of a distorted thermal population represented 
by a Laguerre polynomial expansion. The first two coefficients are A q = 1, Aj 
= 0 with the result that T th = T, and A 2 lies in the range from 0 to .5. The five 
temperatures are shown in Fig. 8 for T = 1000° K. A large temperature difference 
results for values, A, >0.1. Again the three probe temperatures lie above T , 
while the radar temperature lies below. 

An examination of the distribution function in Fig. 9 shows that the slope 
becomes steeper as A 2 increases at low energies with the opposite behavior 
at high energies. A steeper slope results in a smaller temperature. Therefore 
since the radar backscatter samples only the low energy electrons, it gives a 
lower temperature as A 2 increases while the probes sample the higher energy 
electrons and consequently give a higher temperature as A 2 increases. 

The parameter, L, giving the degree of deviation from equilibrium ranges 
from L = .137 for A 2 = .1 to L = .685 for A 2 = .5. Figure 10 shows the varia- 
tion of the relative probe to radar temperature difference versus the extent of 
deviation from equilibrium for the distorted thermal electron population. The 
points from the two-temperature distribution and the polynomial expanded dis- 
tribution tend to follow the same curve. From the curves we see that for 
example a 50% temperature difference corresponds to L values of approximately 
.15 to .34. 

It must be noted that in contrast with the photoelectron distribution function 
we do not have realistic values for the parameters of the distorted thermal 
distribution functions that would typify the ionosphere. Therefore rather than 
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setting limits on what temperature differences can exist in the ionosphere due to 
a known, distortion of the low energy electron distribution we have demonstrated 
what distortion is necessary to produce a given temperature difference. 

SUMMARY 

An examination of the probe and radar electron temperatures for an isotropic 
non-Maxwellian electron distribution function has demonstrated the distinctiveness 
of each probe and radar temperature parameter, in both the energies and the 
methods of sampling the electrons. We find that a realistic photoelection popula- 
tion does not produce significant temperature differences. Large temperature 
differences, 50% or more, can result from higher than normal fluxes of 
energetic electrons or from a distorted low energy electron population. If the 
electron distribution function can be characterized by a different temperature 
above and below about 1KT as illustrated in Fig. 9 then the probe and radar tem- 
peratures will be quite different since they sample in the higher and lower energy 
regions respectively. In the development and use of the probe and radar tech- 
niques, the Maxwellian form of the distribution function has been a basic assump- 
tion. Consequently neither technique has been operated in a manner allowing it 
to measure the electron distribution function at low energies. It is hoped that 
future investigations, both experimental and theoretical, will study the ionospheric 
electron distribution function, particularly in the energy range below 2 ev. Such 
studies will be of great interest not only for the temperature discrepancy problem 
but also for understanding the myriad of ionosphere electron interactions. 

If the isotropic ionospheric electron distribution function were sufficiently 
distorted from a Maxwellian form to produce large probe - backscatter tempera- 
ture differences as illustrated in the sample calculations, then one would not be 


justified in using a single temperature parameter in all calculations of physical 
quantities (such as energy, thermal conductivity, etc.) nor would one be able 
to assign a preference for one instrument temperature over another. The probe 
and backscatter measurements would in this case provide complimentary infor- 
mation. However the relevant parameter to be measured would then be the 
electron distribution function itself. 
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f (x) = e”“ [1 + A 2 {15/8 - 5/2 x + x 2 /2}]. 


Figure 10. Relative temperature difference between probe and radar 
correlated with the extent of deviation of electrons from 
equilibrium. Solid curves computed with polynomial expanded 
dist r '.i-..tion, dashed curves computed with two temperature 
distribution. 




Figure 1 Photoelectron regions for various tem- 
perature and flux values. Numbers at 
70,000 °K represent ambient density in 
particles per cc. 







photo = pyr 2 


Figure 2. Contours of relative difference between 
thermodynamic temperature and for 

photoelectrons. T, = 1500°K. 
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Figure 5. Contours ofrelative difference between 
cylinder probe temperature and T^ for 
photoelectrons. T 1 =1500° K, V =.5 to 
2.5 volts. 



Figure 6. Variation of the-five temperature parameters 
with photoelectron temperature. Tj =1500°K, 
Photo =1. 
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